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Alas  tract 

Tills  study  deals  with  au  electronic  control  system  rodifi cation 
to  solve  a torque  limitation  probl.em  in  i.he  Airiorne  Dynarac  Alignment 
System  aboard  the  Airborne  La.ser  Lab.  Tills  control  system  is  reduced 
in  complexity  and  modeled  in  the  frequonciy  suid  state  space  domains. 
Closed  loop  system  response  over  a desired  frequency  range  is  obtained 
from  the  frequency  domain  model.  Stability  data  is  calculated  through 
eigenvalue  analysis  of  the  state  space  model.  A low  pass  filter  is 
introduced  to  decrease  the  demand  on  the  torque  saturated  motor.  'The 
main  result  if  this  modification  is  that  the  filter  does  not  appre- 
cia,bly  alter  the  closed  loop  response.  This  is  due  to  the  fact  that 
the  cutoff  frequency  of  the  filter  must  be  extremely  high  to  keep  the 
modified  system  stable.  I t is  concluded  that  the  complimentary  filter 
does  not  soive  the  torque  satura.tion  problem.  A recoirimendation  is 
ricde  that  f’urther  analysis  be  concentrated  in  the  time  doniain  through 
modem  state  space  tcctniques  and  the  design  approach  to  an  optimal 


controller  is  presented. 


I.  Introduction 
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Background 

The  backjjroiund  of  the  mirror  alignment  r.yotom  lies  in  the  design 
progression  of  the  Aii'boi'nc  Laser  Lab  (ALL)  based  at  Kirtlauid  AT.B, 
Albuquerque,  New  I'exico.  The  mission  of  the  ALL  is  to  demonstrate 
the  tecimology  involved  in  firing  a laser  from  an  airez'aft  in  an 
actual  flight  environment.  This  project  was  divided  into  several 
phases.  The  initial  phases  involved  the  design  and  fabrication  of  the 
individual  components  aboard  the  ALL.  In  the  final  piiases  all  the 
components  will  be  installed  in  the  aircraft  and  tested.  I mried lately, 
the  need  for  certain  subsystems  to  the  ALL  was  acknowledged. 

These  included  the  actual  team  generation  mochaniSi:i:  spool j i csi iy 
the  Gas  Dynamic  Laser  (GDI),  and  a highly  specified  Autoriiatic  i'cair 
Pointing  and  Tracking  (APT)  device.  Dpon  further  analysis,  however, 
the  need  for  an  active  bean  alignmona  cystei’’.  betwee;,  ohe  Jasei’  and  the 
pointing/tracking  device  became  apparent.  The  main  catalyst  behind 
this  need  was  the  unique  vicrational  environment  encovrtered  in  a 
flying  aircraft. 

The  ALL  is  a highly  modified  XC-133  (Boeing  70?)  aircraft.  An 
aircraft  of  this  size  would  have  significruit  structural  vibrations 
that  would  cause  misalignments  in  the  iorer  .:ean.  This  is  especially 
critical  due  to  the  fact  that  the  ALL  is,  carry! jig  the  GDL  which  is 
primarily  a turbojet  engine.  Tliorefore  the  vibration  of  the  GDL 
coupled  with  i.h.e  ncr'.\ai  twisting  rir.u  bending  of  the  aircraft  present 
quite  an  adverse  vibrational  environment  to  thi  laser  beam.  Thus,  the 
Airborne  Dynamic  Alignment  System  (ADAS)  was  designed  to  eliminate  any 
beam  misalifqiment  due  to  structural  flexures  of  the  ALL  and  the  vibi'a- 
tion  of  the  GDL. 
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Problem 


Hie  main  problem  in  der.igning  the  ADAS  was  quantitatively  defining 
the  vibi-ational  environment.  A considerable  amount  of  information  on 
the  structural  vibration  of  an  ordinaiy  KC-135  is  available.  Nonetheless 
the  ALL  was  not  an  ordinary  airplane.  For  example,  major  structural 
modification  was  necessary  to  install  the  GDI.  and  a hole  was  cut  in 
the  roof  of  the  aircraft  to  mount  the  APT.  Also,  there  was  little 
information  on  just  how  the  GDL  vibrations  would  couple  with  the 
aircraft  vibrations . 

V,'ith  all  tills  in  mind,  a "first  cut"  approximation  of  the 
vibrations  was  derived.  A complex  mathematical  model  of  the  ALL  was 
obtained  and  a "Na.stran"  based  (Ref.  2)  structural  analysis  ’was 
pei’forned.  The  result  of  this  effort  was  the  ADAS  Statement  Of 
Vork  (SCVl)  Power  Spectral  Density  (pSD)  . (Ref.  J-,2.29  and  see  Fig.'s 
1 and  2)  . It  was  with  this  ADAS  SOV/  disturbance  spectrum  that 
Perkin-Elmer  designed  the  actual  ADAS  hardware.  (Ref.  6).  The  system 
was  delivered  to  the  Air  Force  to  be  used  in  the  ALL.  It  was  tested 
and  found  to  be  within  specifications.  Using  a small  laser  for 
testing,  the  ADAS  and  APT  wore  Installed  in  the  aircraft.  Flight  tests 
revealed  that  the  actual  vibrational  en/ironment  was  worse  than  pre- 
dicted. There  were  significant  large  magnitude,  high  frequency 
vibrations  that  were  not  present  in  the  SOV/  PSD  for  ADAS.  Further- 
more, those  large  m2ignitude  high  frequency  disturbances  were  causing 
ADAS  to  saturate.  Essentially,  the  system  could  not  respond  fast 
enough  to  null  the  beam  misalignment.  Since  this  knowledge  was  a 
result  of  CYCLE  II  testing,  the  new  vibration  PSD  became  the  CYCLE 
III  PSD  (see  Fig.  3) • 
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and  a decision  was  made  to  improve  ADAS  to  meet  this  new  input 
disturbance  PSD. 


General  Problem  Solution  Approa.ch 

I 

[ The  reason  the  ADAS  saturates  in  the  presence  of  the  high 

frequency  disturbances  is  the  torque  limit  of  the  mirror  actuator  motors. 
In  this  report  it  is  proposed  to  accept  the  torque  limit  that  exists 
and  correct  the  problem  through  electronic  control  system  modification. 

i Tliis  is  possible  because  the  present  translation  error  is  very  low  and 

more  error  may  be  tolerated.  The  general  approach  to  find  a solution 
is  as  follows: 

f 

1 1.)  Conduct  an  indepth  analysis  of  the  pres- 

I ent  ADAS. 

1 I 

! 2.)  Quantitatively  define  the  extent  of  the 

1 saturation  problem. 

K 3*)  Mathematically  simplify  the  pj-esent  sys- 

I tern  for  ease  of  analysis  and.  design. 

4. )  Modify  the  control  system. 

5. )  Examine  results  of  above  modification. 

Finally,  conclusions  will  be  dravm  from  the  above  investigation  and 
recommendations  for  further  action  will  be  made  based  on  those  con- 

I elusions . 

Limits  of  Study 

The  limits  of  this  study  ar'e  as  follows: 

1. )  Modification  will  deal  only  with  the 

electronic  control  system,  and  no  hard- 
ware changes  were  investigated. 

2. )  Frequency  domain  analysis  will  be  limited 
to  using  the  ADAS  frequency  domain  simu- 
lation . 
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3. )  All  electronic  non-linearities  in  the 

frequency  domain  simulation  of  ADA3  are 
ignored. 

4. )  No  attempt  was  made  to  examine  the  inter- 

active effects  of  ADAS  modification  on 
the  APT  or  the  Auto  Alignment  System 
(aas) . 

Ihe  non-linearities  in  ADAS  are  in  the  form  of  a current  limit  in 
the  compensation  networks  and  a torque  limit  to  each  mirror-acctuator 
motor.  Ignoring  the  limits  will  allow  use  of  a completely  linear 
representation  of  the  control  system.  The  RMS  value  of  the  signal 
in  the  compensation  networks  and  to  each  m.otor  will  be  computed 
during  the  simulation.  In  this  viay,  violations  of  the  limits  and  the 
resulting  saturation  of  the  control  system  will  be  identified. 
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II. 


Analysj  s of  Prei-xm-t,  ADAS 


General  Description 

The  Airborne  ^namlc  Alignncnt  S^ystem  (ADAS)  is  a subsystem  of 
the  Airborne  Laser  Lab  (ALL)  and  serves  to  maintain  alignment  between 
the  Gas  ^namic  Lasei'  (GDL)  and  the  entrance  of  the  Automatic  Pointing 
and  P'acking  device  (APT).  The  ADAS  is  mounted  to  the  .APT.  In  this 
way  there  is  no  misalignment  problems  betv/een  the  API  and  ADAS.  Tne 
two  sources  of  beam  vibration  or  "jitter"  come  from  the  GDL  and  the 
flexure  of  the  ALL.  Since  the  ADAS  is  not  hard  mounted  to  the  GDL, 
the  aircraft  flexure  is  considered  the  lar'ger  bean  misalignment  source. 

The  ADAS  included  the  following; 

A.  tran.slaLion  sensor  (TS) 

B.  prealignment  sensor  (P.AS) 

C.  beam  angie  sensor  (|BAS) 

D.  two  bean  steering  mirrors  (B'Di) 

E.  associated  power  souj'ces  ouid  control  electrcnics 

H)e  PAS  and  Tii  each  have  their  own  active  light  source.  (Ref.  4). 
Tlais  is  to  allow  the  use  of  the  PAS  and  T3  before  the  act\ial  laser 
beam  is  generated.  Tne  BAS  is  used  only  during  the  operation  of  the 
GDL.  It  maintains  a constant  "sa.mple"  of  the  precise  angular  position 
of  the  bean.  All  three  sen.'-.ors,  PAS,  TG,  and  BAS  are  utilized  to  sense 
misalignments  from  a nominal  posltio!i  ajid  feed  this  information  to  the 
electronic  control  system.  The  translation  sensor  measures  transla- 
tional misalignment  in  the  ADAS  optical  train  between  the  GDL  and  the 
traJislation  sensor  Fold  Kd7'ror  (FM)  . The  prealignment  sensor  mea.sures 
angular  alignment  between  the  GDL  and  its  fold  mirror  (PAS  FM) . Ihe 
beam  angle  sen.nor  gets  its  sample  of  the  main  beam  from  the  High 
Power  Grating  (HPG)  and  senses  the  precise  angle  of  the  beam 
between  the  grating  and  the  GDL. 
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Fig;.  4 . Functional  Dia^rran  of  the  C-ptical  Train  and 
Components  of  the  Airborne  Dynamic  Alignment 
Assembly 
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The  electronic  control  system  uses  this  sensor  information  to 
command  the  movement  of  the  two  beam  steering  mirrors.  The  two  beam 
steering  mirrors  are  mounted  to  an  optical  bench  which  is  in  turn 
mounted  to  the  floor  of  the  ALL.  Tije  stop  to  stop  angular  displace- 
ment of  each  miiror  is  10  nicroradians.  Vibration  isolation  mounts 
are  used  on  the  bench  and  GDL  to  effectively  reduce  their  vibrational 
coupling.  The  functional  layout  of  the  ADAS  is  shown  in  Fig.  4. 

The  Electronic  Control  Sysrem 

The  control  system  works  in  essentially  two  independent  cha.nnels. 
Those  arc  the  Y-ALPHA  chaiuiel  and  the  Z-BETA  channel.  Y and  Z are  the 
two  translational  ei-rors  and  ALPHA  and  BETA  are  the  two  angular  errors. 
Each  of  the  two  mirrors,  has  two  sejiarate  axes  of  rotation.  Two 
sepau'ate  mirrors  are  necessary  because  it  is.  not  phyrdcally  possible 
to  null  out  angular  error  without  inducing  tiauislational  err:r  with 
only  one  mirror.  (Ref.  3:2.4,  and  Appendix  B)  . Hie  block  diagra.m  of 
the  control  system  is  shoiai  in  Fig.  5*  Each  of  the  four  loops  is 
functionally  the  same.  A representative  single  loop  is  shoim  in  Fig.  6. 
Each  loop  basically  consists  of  a compensation  network,  electronic 
decoupling,  mirror-motor  dynamics,  optical  cross  coupling,  arid  unity 
feedback , 

The  electronic  cross  coupling  is  necessary  to  decouple  the 
physical  (optical)  cross  coupling.  Essentially,  it  is  with  this  elec- 
tronic decoupling  networks  that  translat-ional  error  is  not  induced 
while  nulling  a pure  angular  disturbance. 
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Reductj  on  of  Control 


stem 
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Reduction  of  the  four  loop  control  system  to  a two  loop  system 
is  straightforward  considering  the  sensitivity  coefficients  between 
the  control  loops  (Fig.  5) • Furthermore,  this  will  reduce  the  difficul- 
ty in  understanding  and  handling  the  control  system. 


Table  I 

Comparison  of  Ei-ror  Output  Between  4 Loop  £ind  2 Loop  Control  System 

4 Loop  System  2 Loop  System 

Y-TRANS  Error  4.l6^^  inches  4.165  inches 

ALPHA-P  Error  1.707  radians  1.705  radians 


In  Table  I,  the  results  of  the  reduction  are  shovfn.  This  is  the 
error  output  with  the  ADAS  SOV/  input  at  the  ALPHA-?  disturbance  point 
with  no  modification  to  the  control  system.  The  tvjc  loop  sysiefi 
gives  almost  exactly  what  the  four  loop  system  gives.  Therefore, 
the  two  loop  system  (Fig.  7)  will  be  used  in  place  of  the  four  loop 
system.  Essentially,  the  Y-TRAH3  — ALPHA.-?  loops  are  exactly  Lhe 
same  as  the  Z-TOANS  — EETA-P  loops. 

Need  for  Electronic  Decoupling 

nie  fact  that  there  is  optical  coupling  between  the  two  beam 
steering  mirrors  necessitated  electronic  decoupling  between  the  con- 
trol loops.  Each  mirror  has  an  ALPHA  and  a FETA  axis  of  rotation. 
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They  are  defined  ro  Uiat  ALPHA  of  mirror  "1  couples  wilh  ALPHA  of  ! 

1 

mirror  H2  and  the  nai:'.e  for  the  FETA  axis.  The  electronic  decoupling 
network  for-  the  Y-TRANLLAnCH  and  the  ALraA-FHEALIGN  control  loops  is 
shown  in  Fig.  8.  Tins  cross-coupling  is  of  the  sane  form  as  thar 

between  the  Z-TRAHSLATIOM  and  EETA-FREALIGN  control  loops. 

' 

i 


KYT  KAIYT 


i 


rg  • 8 . .vsi  liA— i f 


i'..AL.j  ..Icc ironic  Grosscoupling  Hetwork 


The  Y-TRANS  and  ALPHA-P  control  loops  are  shown  septvrated  from  the 
main  clock  diagram  in  Fig.  9-  In  thi:-.  figure,  however,  only  the 
V constant  gains  and  sensitivities  are  included.  U.e  electronic  cross- 

coupling sensitivities  (KAl.A?  and  KA2YT)  are  labeled  A and  B respec- 
tively. The  high  frequency  blocks  are  also  not  Included  for  further 
simplification.  In  this  way  the  control  loops  are  as  they  would  be 
at  a very  low  frequency  (s-^O)  condition. 

The  equations  for  each  error  will  be  written  and  the  values  of  A 
and  R will  be  found  such  that  a pure  angular  disturbance  (a(I)IS))  will 


I 
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not  induce  a trannlatiotiai  disturbance  (T(Dirj))  . Equations  deecrib- 
ing  the  errors  E(t)  and  E(a)  ai’e: 

E(T)  = T(DIS)  +2.41  ALPHA(l)  - 5-05  ALFKA(2) 

E(a)  ---  A(DIS)  + 3.998  ALFHA(2)  - 3.86  ALrKA(l) 

The  equations  describing  ALFHA(l)  and  ALrHA(2)  are: 

ALFHA(l)  = (12)  (7.1E-4)  (2.78)  (3000)  (E(t)) 

+ (B)  (2.78)  (5000)  E(a) 

ALFHA(2)  - (12)  (7.1E-4)  (0.477)  (15.3)  (5000)  (E(A)) 

+ (B)  (2.78)  (5OOG)  E(t) 

Substituting  ALFr.A(l)  aj'.d  ALrhA(2)  ivitc  the  E(T)  aquation; 

E(T)  =:  T(DIS)  + (2.41)  (118.4)  E(t)  + A(631.S)  E(a) 
-5.05  (310.9)  E(a)  + B(118.4)  E(t) 

= T(DI3)  + (2S5.3)  E(t)  + A(lq70.8)  E(a) 

- e(598)  e(t)  - (1570.05)  e(a) 

= T(DIS)  + (285.3  - 5^SB)  E(T)  + (15?0.8A 

- 1570.05)  e(a) 

Now  in  order  to  make  E(T)  a function  only  of  T(DIS)  and  not  of  E(a): 

(1570.8A  - 1570.05)  must  equal  0 
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Tl'ierefore: 


A = (1570.05)/(1570.8)  = 1.0 

Ey  the  same  process  B is  found  to  he! 

B = 0.965 

In  the  original  block  diagram  (Fig.  5))  A = 1.0  and  B = 0.95»  so  the 
electronic  cross-coupling  (decoupling)  does  Indeed  decouple  the  angle 
and  translation  control  loops  at  low  frequencies. 

Discussion  of  Error  Limits 

The  error  limits  specified  in  the  ADAS  SOW  are: 

Table  II 

Error  Limits  In  Angle  And  Translation 
Signal  Liml t 

ALFHA-P,  BETA-P  3-0  radians 

Y-TRANS,  Z-TOALS  2,500.0  inches 

In  the  interest  of  keeping  this  report  unclassified,  the  exact 
scale  of  the  limits  and  cloned  loop  response  is  not  exact. 


L 


18 


Ill . ADAS  Frequency  Doir.aln  Simulation 


ADAS  ip  modeled  in  the  fx'equency  domain  throuch  a general  purjicse 
program  called  FRQRSP.  FRQHSP  was  developed  at  Oklahoma  State  Univer- 
sity (Ref.  4)  for  the  Air  Force  Weapons  Laboratory.  FRQRSP  solves  a 
system  of  simultaneous  linear  equations  by  the  Gai'ss- Jordan  elimina- 
tion method  over  a specified  frequency  range.  The  program  is  capable 
of  accepting  input  info3.'mation  in  several  modes  including  distui'bance 
spectrum  data  Siich  as  the  ADAS  SCR  and  CYCLE  III  PSD. 

ADAS  was  modeled  by  a set  of  l6  equations,  each  representing  a 
particular  node  in  the  control  system.  The  l6  equations  are  written  :j 
the  matrix  form: 


11 

- 

1 

A 

H 

B 

J 

L 

where 


[■]• 

[■]■ 


is  the  vector  of  the  unknown  states 

is  the  vector  of  inputs  to  the  system 
is  the  matrix  that  describes  the  system 


As  aji  oxanple  of  the  way  in  which  the  system  is  modeled,  examine 
the  equation  for  Y-Translation  error.  From,  Fig.  7: 


Y-TRAN3  Error  = Y-TRANS  Disturbance  - (-2.4l)  ALPHA(l) 

- (.5.07)  ALFHA(2) 
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where: 


ALPHA(I)  is  the  An{;le  of  Mi-^rroi  til 
ALPHA(2)  is  the  Angle  of  Hirror  #2 

Rearranging  the  equation: 

Y-TRANS  Error  + (-2.4l)  ALPHA(I)  + (5-07)  ALPKA(2)  = Y-TRANS  Disturbance 


In  matrix  form: 


' 

- 

1 -2.41  5.07 

Y-TRANS  Error 

- 

1X3 

ALPKA(l) 

Y-TRANS 

ALPHA(2) 

- 

Disturbance 

- 

3X1 

. 

Thus  the  Y-TRAIIS  Error  equation  is  in  Ibie  form: 


- — 

- 

A 

R 

B 

.. 

As  an  example  in  which  the  frequency  sensitive  portions  of  the 
control  system  are  modeled,  examine  the  compensation  network  in  the 
Y -Translation  control  loop: 
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Let  the  signal  at  Node  it2  equal  N0DE2SIG  and  let  the  Y-Tx'anslation 
error  signal  equal  YTERSIG.  Thun: 


N0DE2SIG 


(YTERSIG)  (5000) 


9s 

3140  +"0 


f-  1.0 


= YIERSIG 


45,000s 
3l40  + s 


+ 5000 


= YTERSIG 


■ 1.57E  + ? + (5.OE  + 4)s 
3l40  + s 


Rearranging  this  equations 


1.57E  + 7 + (5.OE  + 4)s  ■ 
3l^i0  + s 


- 

— — 

YTCRSIG 

N0DE2STG 

— 

— « 

which  is  in  ihe  correct  form  of: 


- ^ 

r 1 

A 

B 1 , 

B 

..  ^ 

L J 

This  procedure  is  followed  for  the  other  14  nodes  and  then  ail  the 
equations  ai’o  combined.  The  result  is  a system  of  equations  in  the 
form: 


“1 

■ 

* 

A 

R 

= 

B 

J 

»■  •. 

16X1 

16X1 


Since  the  frequency  response  of  the  system  is  desired,  replace  s 
by  jw.  Ttien  evaluate  at  a paxticular  frequency  w and  the 

result  is  a complex  [A],  The  complex  [tl  matrix  is  then  inverted 
and  multiplied  vrlth  to  result  in  a complex  [^rJ  . The 
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procedure  is  then  repeated  for  a certain  range  of  frequencies.  The 
overall  result  is  then  the  complete  system  response  over  that  certain 
reinge  of  frequencies. 

The  computer  simulation  model  of  ADAS  is  limited  in  that  it  is 
not  possitle  to  model  the  nonlinearities  such  as  motor  torque  limits 
and  current  limits.  It  is  also  not  possible  to  program  multiple 
inputs.  Hovfovor,  once  the  system  of  equations  is  solved,  output  may 
come  from  any  node. 
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IV . ADAS  State  wpacc  Model 


ADAS  is  modeled  in  state  space  form  to  satisfy  the  need  for  more 
precise  stability  information.  Itiis  section  will  develop  the  state 
equations,  and  the  shaping  filters  designed  to  simulate  the  input 
disturbance  P3Ds.  Tl'.e  states  of  the  control  system  are  labeled  in 
Fig.  10.  In  this  particular  figure,  only  the  Y-TR/JIS  and  ALPHA-P 
control  loops  are  shovm.  The  figure  would  be  exactly  the  same  for  the 
Z-TRANS  ESTA-P  control  loops. 

The  block  diagram  for  each  motor  is: 


1^10? 

140?  + s^ 

--  Y 


The  state  equations  for  each  motor  is: 


where: 


— 

~ — 

0 1 

0 

X + 

-140?  0 

1407 

— ~ 

, for  motor  01 


, for  motor  02 


Kotor  #1 


Fig.  10.  Location  of  States:  in  ADAS  Two  Loop  Control  Systen 


The  actual  block  dlacran  for  the  Y-TR/u\’S  (GGYT)  compensation  io: 


I 


Tlie  X, 
j 


state  equation  is; 


and 


The  input  disturl'^^nces  must  be  approximated  by  white  noise  passed 
through  a shaping  filter.  The  shaping  filter  for  the  ADAS  SOW  PSD 


I 
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for  the  Y-TRJi'_^’3  control  loop  is; 


where 

W = white  noise,  zero  mean,  Gaussian  statistic 

K = O.l^iE-10 


K 


+ Bs  + C 


B = 7.5+E+2 
C = 1.1I2E+2 


The  shaping  filter  for  the  ADAS  SCW  PSD  for  the  ALPKA-P  control  loop 
disturbance  is: 


where 


W = white  noise,  zero  mean,  Gaussian  statistics 

K = O.llOE-8 


Xq  K 

✓ 

_ ^ O " 

w s + Bs  -f  Cs  + Ds  + S 

B = l.OE+4 

C = 3,78E+5 

D = 6.34E+7 

E = 3.97E-(9 

iTie  state  equations  are; 


0 

1 

0 

0 

1 

o 

L 

0 

0 

1 

0 

0 

0 

0 

0 

1 

X + 

0 

-E 

-D 

-c 

-B 

K 

— 

— ^ 

where 
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f 

I 

7 

f 

I 


L 


The  origin  of  these  shaping  filters  is  hotter  understood  by  plotting 
the  riagnitude  versus  frequency  curves  for  the  ADAS  SOW  and  approximte 
disturbances.  These  plots  for  the  Y-TRAXS  and  ALPHA-P  control  loops 
are  shown  in  Fig.’s  2?  and  28  respectively. 

The  coLihined  state  space  ^aJ  natrix  is  shovni  in  Pig.  24  in 
Appendix  C.  The  eigenvalues  of  this  |^Aj  matrix  coz*respond  to  the 
roots  of  the  characteristic  equation  of  the  closed  loop  transfer 
function.  The  eigenvalues  for  this  control  system  are  shown  in 
Table  III. 


Table  III 

Eigenvalues  of  AD.^.S  Control  System 
(2  Loop) 

-4184  + 844i 

-1497  + lOSfci 

-897 

-179 

-388 

-365 

-202  + 24i 
-297  + 83i 

The  fact  that  all  eigenvalues  are  negative  means  the  original  ADAS  is 
stable.  (Ref,  5: 576) . A listing  of  program  EIGEN  that  solves  for  the 
eigenvalues  and  eigenvectors  of  ^aJ  is  in  Appendix  C.  Consult  Ref.  7 
for  explanation  of  subroutine  RGEIG  used  in  computer  program  EIGEN. 
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V.  Nur.crlcal  Kxtcrit  of  Ihe  FrobJon 


The  baseline  perfornance  for  ADA;'  using  the  Statement  Cf  Work 
(sow)  distui'bance  is  given  in  Table  IV.  Those  basleine  figures  will 
be  used  in  comparsion  to  the  worst  case  disurbance  and  also  the 
results  of  any  compensation . 

As  noted  in  the  introduction,  the  worst  case  disturbance  is 
called  the  CYCLE  III  disturbance.  Since  the  angular  errors  are  of  a 
much  greater  ccnccrn  than  the  translational  distui'bances,  the  CYCLE  III 


Baseline 
(input  is 

Signal 

A1  Torque  Required 
A2  Torque  Required 
Y -TRAN 3 Error 
ALFfiA-P  Error 


Table  IV 

Closed  Loop  Response 
ADAS  SOW  at  Node  -vp) 

Output  (RMS) 
1.29  in -lb 
0.62  in-lb 
4,i6  in 
1.71  rad 


ROLL  RACE  MOTION  PCD  will  be  used.  Due  to  the  physical  layout  of  the 
ADA.3  in  the  ALL,  the  roll  base  motion  aireetly  affects  the  angular 
control  loops.  Similarly,  pitch  base  motion  affects  the  translational 
control  loops.  It  is  for  this  reason  that  for  a single  input  the 
worst  case  system  output  rjcurs  when  CYCLE  III  ROLL  BACK  MOHON  PSD 
is  input  at  the  ALFHA-P  distur’oance  point.  iTie  results  of  this  worst 
case  are  listed  in  Table  V. 
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Table  V 


Worst  Case  Closed  Loop  Response 
(input  is  CYCLE  III  ROLL  PSD  at  Node  //5) 

Output(RHS) 

4.73  in-lb 

2.74  in-lb 
5.30  in 
3.42  rad 

Observe  that  for  both  the  baseline  and  the  worst  case  input,  the 
torque-required-ratio  of  motors  A1  and  A2  is  two  to  one  (2:1).  Under 
the  CYCLE  III  Roll  input  this  two  for  one  torque  demand  exceeds  the 
torque  available  for  motor  A1 . Tliis  condition  occurs  due  to  the 
high  frequency,  large  amplitude  disturbances  in  the  CYCLE  III  ROLL 
PSD. 

Of  equally  great  concern  is  the  increase  in  ALPHA-P  error  over 
the  baseline  figure.  The  ALPHA-P  error  caused  by  the  CYCLE  III  ROLL 
PSD  (worst  case  is  well  over  the  specified  limit,  see  Table  TI  ) . 

Since  these  figures  are  derived  from  a linear  simulation  of  ADAS,  the 
actual  system  would  have  an  even  higher  ALPHA-P  error.  This  is 
because  as  the  torque  limits  are  exceeded  the  control  system  would 
saturate  and  the  motors  would  no  longer  respond  to  correct  for  errors. 

Documentation  is  also  necessary  in  the  form  of  ALPHA-P  open  loop 
response.  The  modification  will  directly  affect  the  open  loop  response 
of  the  angular  control  loop.  Therefore,  it  is  important  to  know  the 
response  of  the  original  system.  This  data  is  shown  in  Table  VI. 


Signal 

A1  Torque  Required 
A2  Torque  Required 
Y-TRAN3  Error 
'ALPilA-P  Error 
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Table  VI 


Open  Loop  Daia  for  ALPHA-P  Control  Loop 

ALPIIA-P 

Gain  Margin  +3.4  dB  3 920  Kz 

Phase  Margin  ■>■2?.  degrees  ij  725  Hk, 

0 dP  Gain  Crossover  725  Hs 

The  magnitude  and  phase  angle  plot  is  shown  in  Fig.  33  • Since  the 
oi)en  loop  staid  liLy  response  is  inaependent  of  the  magni  iudo  of 
disturbance  input,  thoi'e  is  no  need  for  a bascline/worst  case 
compai'ison . 

Ihc  ALPHA-P  open  loop'  Lransfer  function  is  defined  as: 

angular  position  of  laser  beam 

G(n)  = 

error  in  angulac-  alignment  of  the  laser  beam 

The  error  in  angulai’  alignment  is  measured  at  Node  in  Fig.  7 and 
the  angular  position  of  the  bean  is  Node  + /■•^4)  . The  op'on  loop 
roj-.p'onse  is  measured  this  way  to  ensure  unity  feedback  between 
Mode  (/!'8  + tth')  a.nd  the  disturbance/errer  sur.rjng  i)oint. 
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VI . Control  System  Modification  and  Pj'cdlctionr; 

Addition  of  Low  Fass  Filter 

Ihe  control  system  is  modified  ty  piacin-s  a low  pass  filter  in 
the  olectz’onic  cross-coupling.  In  the  frequency  domain,  this  filter 
replaced  the  clock  labeled  KAIAP  (see  Fig.  8)  in  the  oranch  coming 
from  the  ALFrIA-F  control  loop  to  the  Y-TRANS  conti'ol  loop  (Fig.  11)  . 
The  filter  is  named  GLOV.TS,  and  is  of  the  form: 

GLOWPS  = 

(1/T)s  + 1.0 


GLOViTS  is  essentially  a filter  with  T as  tlie  cutoff  frequency. 

A first  order  filter  was  chosen  because  of  simplicity  and  the  desire 
not  to  introduce  more  poles  into  the  system  than  necessary. 

The  cutoff  frequency  T was  varied  between  0 and  700  Hz.  Closed 
looi>  response  data  was  calculated  at  each  frequency  to  compare  errors 
and  required  torques.  Open  loop  data  for  tlie  ALPHA-P  loop  will  be 
gathered  for  stability  considerations. 

In  the  state  space  model  tlie  complimentary  fii+er  will  add 
another  state  to  the  system,  liie  state  matrix  will  be  of 

dimension  13X13  for  the  modified  control  system.  Again  the  filter  wj.ll 
replace  the  block  labeled  KAIAP.  H'.e  state  space  repiesentatlon  of  the 
filter  is: 


j 
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where 


L 


13 


[-T] 


13 


[T] 


substituting  expression  for  U, 


^13 


-(T)  - (T)  -4.4E 


4E+6  X, 


+ 4.51E+6  X^  + (15.3)  X^  + 1.13E+6 
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The  addition  of  state  X^^  also  changes  the  X^  equation.  The  form  of 
X^  equation  for  the  modified  control  system  is; 


X^  = -1407  X^  + 1407  U 


-1407  X, 


+ (1407) 


(12)  (7.1E-4) 


— 

(2.8C^)  -f  12  X^^ 


4.03E-t'6  X,  - 8.5E+6  X^  + (33.6)  X,  + i.6£'E-H6  X, 

I J J f 


+ 12.0 


The  full  state  matrix  for  the  modified  control  system  is  Fig.  22 

in  Appendix  C. 

The  function  of  this  filter  is  to  pass  the  low  frequency  content 
of  the  signal  crossing  from  the  ALPHA-P  control  loop  to  the  Y-TiiAN3 
control  loop  and  attenuate  the  hifher  frequency  portion.  Thus,  the 
intent  of  the  filter  is  to  cut  down  on  the  torque  requirement  command 
of  Kotor  #1.  The  big  difference  between  the  ADAS  SOVJ  and  the 
CYCIJ:  III  roll  disturbance  is  the  large  ma^itude,  high  frequency  error 
in  CYCLE  III  roll.  The  filter  will  effectively  reduce  the  CYCLE  III 


3^^ 


torque  required  signal  to  something  similar  to  the  ADAS  SO'rf  for 
Motor  id  . 


Reducing  the  Motor  #\  torque  will  destroy  the  2:1  mirror  rotation 
for  a pure  angular  ei'ror.  Therefore,  it  is  expected  that  translation 
error  will  increase  witn  the  addition  of  the  filter.  Since  it  is 
possible  to  correct  fo.r  angular  error  with  one  to  one  novenent  of  the 
mirrors,  the  ajigular  error  is  not  expected  to  increase.  However,  with 
the  physical  and  electrical  cross-coupling  present,  the  direct  efi'ect 
on  the  angular  error  is  not  obvious. 
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VII.  Results 


Tlie  frequency  donain  closed  loop  response  results  with  the 
complimentai'y  filter  in  ADAS  is  shown  in  Table  VII.  The  cutoff 
frequency  is  varied  from  32  Hz  to  700  Hz  and  the  input  is  CYGLil  III 
ROLL  PSD  at  the  ALPHA-P  disturbance  point. 

Table  VII 

Oonparison  of  Conplir.entaiT  f'ilter  Cutoff  Frequencies 
(Output  in  RMS;  1-SIGlLA) 


Cutoff 

Free coney 

A1  Torque 
Reaunred 

A2  Torque 
Rc'Quirod 

AIPHA-P 

Y' TRANS 
Eri’o?: 

(in-lb) 

(in- lb) 

(rad) 

(rad) 

32 

1.90 

2.27 

1.75 

6.20 

56 

1.71 

2.26 

1.90 

6.42 

64 

1.65 

2.25 

1.95 

6.49 

111 

1.25 

2.22 

2.33 

6.87 

143 

1.00 

2.20 

2.64 

7.09 

150 

0.96 

2.19 

2.71 

7.13 

190 

0.88 

2.18 

3.14 

7.34 

200 

0.92 

2.18 

3.25 

7.39 

35C 

2.93 

2.49 

4.90 

7.56 

700 

14.80 

8.49 

11.30 

8.90 
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■nils  frequency  domain  cloned  loop  data  vernus  cuLoff  frequency  is 
graphed  in  Fig.'s  12,  13,  1^,  and  15*  Cn  each  of  the  four  graphs  the 
original  system  value  is  plotted  as  a dotted  line  for  comparison. 

From  the  plot  of  Y-TRAMS  versus  cutoff  frequency,  it  is  significant 
that  for  every  frequency  the  modified  system  yields  a highci*  value  of 
error  than  the  original  system.  This  is  as  predicted.  It  was  not 
predicted  that  the  Y-TRANS  error  would  increase  with  increasing  cutoff 
frequency.  Ihis  result  did  not  seem  likely  due  to  the  modification 
involved.  Ac  the  cutoff  frequency  gets  larger  in  a low  pass  filter, 
one  would  expect  closer  agi'eement  with  the  original  system. 

This  phenomena  may  also  be  observed  in  the  ALPHA-?  error  graph. 
The  low  cutoff  freque:icies  give  the  desired  reduction  in  error  (from 
the  original  system),  but  high  ciitoff  frequencies  give  ihe  reverse 
effect.  From  these  two  plots,  one  may  see  that  the  lowest  cutoff 
frequency  yields  the  least  error.  In  fact,  0 Ks  cutoff 
(no  ALPHA-P — ► Y-lnAN'S  electronic  coupling)  would  appear  to  be  the 
best.  Completely  eliminating  the  Al.PJiA-P  — > "''-TRAk'd  elccti’onic  cross- 
coupling was  not  in  the  original  plan,  tut  the  results  are  shown  ir. 
Table  VIII. 


Table  VIII 


Closed  Loop  Response  With  No  ALPHA-P  — ►Y-lTt\!!d  Electronic  Cross-coupling 


A1  Torque 

A2  Torque 

ALPHA-P 

Y-TRANS 

Required 

Required 

Error 

Error 

2.12  in-lb 

2.28  in-lb 

1 .65  rad 

5.86  in 
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Conplinontary  P'iltor  Cutoff  Froq 


Cutoff 


.Tojj-  a, -A 


onpllnentary  Fjitor  Gutoi 


The  A1  and  A2  motor  torqre  plots  each  yield  a riniriium  at  about 
200  Hz.  The  A1  torque  was  considerably  i-educsd  as  predicted.  It  vras 
not  expected  that  the  A1  motor  torque  requirement  would  fall  below 
that  of  the  A2  motor. 

The  AIjFHA--P  open  loop  results  arc  shovni  in  Table  IX  on  the  next 
page.  This  table  shows  the  increase  ir.  frequency  response  stability 
parcumeters  from  the  original  system  to  the  system  with  no 
ALPHA-P  — .-Y-TRANS  electronic  cross-coupling.  Ti;is  result  reinforces 
the  data  from  the  ALPHA-P  arid  Y-T'RANS  error  plots.  The  190  -*  200  Hz 
cutoff,  which  was  the  minimum  region  for  the  torque  required  curves 
yields  an  unstable  open  loop  phase  margin.  ALPHA-P  open  loop  error 
rejection  curves  for  various  cutoff  frequencies  shown  in  Appendix  C. 

The  state  space  ruialysis  results  in  a set  of  eigenvalues  and 
associated  eigenvectors.  TTic  system  is  considered  stable  if  all  the 
eigenvalues  ai-e  negative.  Tiie  results  cf  the  eigenvalue  analysis  are 
as  follovjs; 

1. )  The  0 Hz  cutoff  frequency  (,no  ALPH.A-P — »-Y-TH.AHS 

electronic  coupling)  proauces  a positive  eigen- 
value. 

2. )  It  is  not  until  the  cutoff  frequency  is  greater 

than  725  Hz  that  all  eigenvalues  become  negative. 

3. )  The  system  does  not  appeal'  to  be  well  damped 

(Damping  Ratio  0 .70?)  until  the  cutoff  frequency 
is  moved  past  5i000  Hz. 

4. )  At  approximately  l6,000  Hz  cutoff  frequency,  the 

modified  system  eigenvalues  cire  within  15/^  of 
the  original  system  eigenvalues. 

The  exact  numerical  results  are  presented  in  Appendix  G,  Table  .XI. 

As  a result  of  the  above  eigenvalue  information  more  frequency 
domain  information  was  needed.  This  time  the  cutoff  frequency  was 


Table  IX 


Comparison  of  ,AL?HA-P  Open  Loop  Characteristics 
with  Various  Filter  Cutoff  Frequencies 


Cutoff 

Frequency 

(Hz) 

0 dB  Gain 
Crossover 
(Hz) 

Phase 

Margin 

(degrees) 

Gain 

Margin 

(dP) 

0 

397 

445 

4-12 

@ 918  Hz 

32 

397 

443 

+9 

@ 775  Hz 

56 

380 

4-40 

+9 

@ 775  Hz 

400 

440 

49 

@ 775  Hz 

111 

140 

4-150 

+7.5 

@ 725  Hz 

150 

460 

-380 

* 

200 

.140 

-351 

* 

330 

550 

4-5 

4-2.1 

@ 600  Hz 

350 

550 

+3 

4-1  .9 

@ 600  Hz 

700 

652 

0 

4-2.5 

@ 623  Hz 

Worst  Case 

725 

4-22 

4-3.4 

@ 920  Hz 

* Phase  angle  never  equals  -180  degrees  below  10,000  Hz 
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* varied  from  800  to  l6,000  The  resultr.  are  fiiown  in  Table  X. 

More  information  concerning  the  erratic  nature  of  the  open  loop  i 

' response  is  shown  in  Fig. 's  29  - 32.  i 

[ ' 


VIII.  Concl  jr.i  on.*^ 

Throe  typoa  cf  tacic  analysis  efforis  were  accomplished  to 
analyze  the  performance  of  the  proposed  low  pass  filter  compensator: 

1. )  Op'on  loo]>  frequency  domain  analysis 

2. )  Closed  loop  frequency  domain  analysis 

3. )  State  space  eigenvalue  analysis 

Thie  conclusions  of  each  approach  are  discussed  separately  and  then 
the  sunmarj'  of  the  overall  conclusions  is  presented. 

Open  Loop  Frequency  Domain  Analysis 

The  ALPHA-P  open  loop  transfer  function  is  defined  as  foll.ows: 

angular  position  of  laser  beam 

G(s)  = 

error  in  angular  alignment  of  laser  beam 

Tnis  is  shown  mere  clearly  in  Fig.  7 . llie  error  in  ang’clar  alignment 

is  Mode  #5  and  the  angular  position  of  the  beam  Nodes  (,¥8  4 tfh.)  . The 

ALPHA-P  open  loop  transfer  function  contains  the  ALPHA -P  compensation 

and  the  dynamics  from  nirror-niotors  and  n2.  The  transfer  function  is 

defined  this  way  to  include  the  optical  coupling  between  the  two  beam 

steering  mirrors  (see  Appendix  B)  and  the  electronic  decoupling  from 

ALPHA-P  to  Y-TRANS  loop.  Since  the  open  loop  response  is  independent  of 

input  the  modified  system  is  compared  to  the  original  system  independent 

of  input  angular  disturbanne.  The  open  loop  results  of  the  original 

system  are  shown  in  Table  VI  and  Fig.  33.  This  is  the  response  when  the 

block  labeled  K/ilAP  in  Fig.  8 contain  a constant  gain  equal  to  1.0. 

To  implement  the  low  pass  filter,  KAlAP  was  reduced  with  1 .0 , 

(iTtFs  4 1 
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where  T is  tiie  cutoff  frequency.  The  cutoff  frequency  j.s  varied  I’rom 
0-*  16,000  Hz  and  the  2’cr.ults  are  shown  in  Table  IV,  Table  X and  Fig.'s 

29-32. 

The  conclusion  is  that  the  modified  system  ini'roved  the  ALPHA-? 
open  loop  stability  with  cutoff  frequencies  below  100  Hz.  '.tith 
cutoff  frequency  above  100  Hz,  the  modified  system  actually  decreased 
the  stability  of  the  control  system.  The  gain  and  phase  mao'gins 
became  very  erratic  as  cutoff  frequency  increased.  Although  it  is  not 
clear-,  the  fact  that  the  ALPHA-?  and  Y-TRANH  loops  are  no  longer 
exactly  decoupled  may  explain  the  instability.  Ttie  poor  stability  data 
could  be  a result  of  the  way  the  ALPHA-?  open  loop  transfer  function  is 
defined.  With  the  optical  coupling  and  electronic  decoupling  between 
the  two  loops,  straightforward  freq-rency  dor -tin  stability  analysis  is 
questionable  if  not  impoc.sibio.  For  this  reason,  the  accuracy  of  the 
improved  stability  figures  at  low  cutoff  frequencies  must  not  be 
trusted.  Clearly  the  stability  analysis  of  the  modified  control  system 
should  not  be  done  in  the  freqiiency  domain. 

Closed  Loop  Frequenev  Domain  Analysis 

Using  the  two  loop  ADAS  control  system  shewn  3n  Fig.  7,  the 
CYCLE  III  Roll  angular  disturbance  PSD  was  input  to  the  ALPHA-?  control 
loop.  The  translational  disturlance  vras  kept  equal  to  zero.  Ihe  four 
main  points  for  measuring  the  closed  loop  response  were:  tlie  signal 
(torque  command)  to  mirror-motors  and  tIZ,  and,  the  Y-TRANS  and 
ALPHA-?  error  signals.  The  results  of  this  simulation  are  shown  in 
Table  V and  are  called  worst  case  figures.  A first  order  low  pass 
filter  is  substituted  for  block  labeled  KAIA?  (see  Fig.  8)  and  the 
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cutoff  frequency  war.  vai'ied  from  0— ^l6,000  Hz.  Tb.e  cloEod  loop 


results  ai'o  presented  in  Tables  VII,  VIII,  ajid  X,  aiid  Fjg.'s  12-15. 

At  low  cutoff  frequencies  (below  200  Hz)  the  modified  system 
closed  loop  res.ponse  was  gi'eatly  improved  over  the  worst  case.  Hie 
torque  requirements  went  down  and  the  ALPHA-P  error  was  reduced.  Hie 
Y-TRAiIS  error  increased  slightly  but  is  still  well  within  the  limits  of 
the  system  (sec  Table  II  ) , Hie  torque  requirements  to  the  Beam 
Steering  Hirror  (BSfl)  //I  and  "2  were  minimized  at  about  200  Hz  cutoff 
frequency.  Howevei',  overall  system  response  seemed  best  at  very  low  or 
zero  cutoff  fx'equcncio.' . The  errors  were  at  their  lowest  value  with 
a zero  Hz  cutoff  and  the  torque  i-equircmcnts  to  each  BSt5  were  well 
below  theij-  worst  case  values.  Using  a 0 Hz  cutoff  frequency  resulted 
in  an  almost  one-for-otic  mirror  torque  rcquii'cmcnt.  Hiis  means  tnat 
both  mirror-motors  .are  Icing  driven  by  exactly  the  same  amount.  There  is 
no  longer  a two-for-one  torque  requirement  between  B3M  ,fl  and  BdH  it'e. 

Hiis  two-for-one  torque  requirement  is  undesirable  because  EUI'I  ffl. 
reaches  its  toi'que  lirdt  and  the  conti'ol  system  saturates.  Table  VI J 
shows  that  the  two-for-one  torque  requirement  is  reduced  to  one-to-one 
at  the  expense  of  increasing  the  Y-HlANS  error.  For  0 Hz  cutoff,  the 
Y-TRANS  error  is  still  vfell  within  the  limit  (see  Table  VIII  and 
Table  II  ). 

However,  as  cutoff  frequency  is  increased  between  250  and  700  Hz, 
the  ALPHA-P  error  exceeds  the  worst  case  error  (see  Fig.  l4) . Since  the 
effect  of  the  low  pass  filter  decreases  as  the  cutoff  frequency  increases, 
the  result  that  the  ALPHA-P  error  exceeded  the  worst  case  was  unantic- 


response  does  approach  the  worst  case  figures.  The  problem  is  that  tlie 
closed  loop  errors  and  torque  approach  the  worst  case  figures  from  the 
high  side.  This  mecuis  that,  above  250  Hz  cutoff  frequency,  the  low 
pass  filter  does  not  provide  any  impz’ovement  in  perfoi-mance.  Also, 
above  about  700  Hz  the  i;;cdifiod  system  results  are  academic  because 
the  system  response  is  not  this  high. 

The  conclusion  from  this  closed  loop  fi’equoncy  domain  analysis  is 
that  the  cutoff  frequency  must  be  below  250  Hz  for  the  low  pass  filter 
to  improve  the  performance  of  the  system.  In  fact,  from  an  error 
control  standpoint,  0 Hz  cutoff  frequency  seems  to  be  the  best.  A low 
pass  filter  vdth  0 Hz  cutoff  frequency  is  the  same  as  completely 
eliranating  the  electronic  cross-coupling  from  the  ALPHA-P  to  the 
Y-TRAMS  control  loops.  This  conclusion  sup'ports  the  belief  that  the 
exact  decoupling  must  be  reduced  or  eliminated  to  reduce  the  torque 
saturation  problem  on  BSH  rl  and  also  control  errors. 

State  Space  Sigenva’iuo  Analysis 

Tne  closed  loop  eigenvalues  for  the  original  and  modified  system.s 
were  calculated.  The  original  system  was  chaincterized  by  12  states 
(see  ADAS  State  Space  Model)  and  tiie  addition  of  the  first  order  low 
pass  filter  added  ai.  additior.al  state.  The  eigenvalues  wore  calculated 
from  the  closed  loop  system  matrix  using  program  EIGEN  (listed  in 

Appendix  C)  and  subroutine  RGEIG  (Ref.  7). 

Analysis  of  the  original  closed  loop  system  yielded  12  eigenvalues 
with  negative  real  parts  (see  Table  III).  This  means  that  tne  original 
closed  loop  system  is  stable.  (Ref.  5: 5?6) . The  eigenvalues  were  then 
calculated  for  the  modified  control  oysto.m  with  the  low  pass  filter 
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cutoff  frequency  vaxled  from  0— ►l60,000  Hz.  The  first  7 eigenvalues 
for  each  case  are  shown  in  Table  XI.  The  last  6 eigenvalues  corres- 
pond to  the  6 states  of  the  disturbance  shaping  filters  (see  ADAS 
State  Space  Model) . These  6 eigenvalues  have  no  influence  on  the 
stability  of  the  system  so  they  are  not  shown. 

The  eigenvalues  of  the  modified  system  show  that  it  is  unstable 
for  cutoff  frequencies  below  about  725  Hz.  As  the  cutoff  frequency  is 
increased  above  725  Hz  the  modified  system  eigenvalues  approach  the 
original  system  eigenvalues.  The  only  difference  is  the  presence  of  an 
eigenvalue  corresponding  to  the  low  pass  filter.  As  cutoff  frequency 
is  increased,  the  effect  of  this  eigenvalue  of  the  filter  decreases. 
This  result  is  as  expected  and  corresponds  to  the  declining  effect  of 
the  low  pass  filter  as  cutoff  frequency  goes  to  infinity. 

Since  the  open  loop  response  of  the  control  system  is  not  greater 
than  725  Hz,  any  filter  with  cutoff  frequency  above  725  Hz  would  have 
no  effect  on  the  system.  And  since  the  modified  system  is  always 
unstable  at  cutoff  frequencies  below  725  Hz,  the  low  pass  filter  is 
not  successful  in  improving  the  system. 

Overall  Conclusions 

Combining  the  conclusions  of  the  three  analysis  techniques,  the 
overall  conclusion  is  that  the  modified  ADAS  control  system  offers  no 
improvement  in  performance  over  the  original  system.  However,  it 
is  interesting  to  compare  and  contrast  the  results  and  conclusions  of 
the  frequency  domain  techniques  with  those  of  the  state  space  analysis. 
Both  the  open  and  closed  loop  frequency  domain  results  pointed  to  an 
improved  stable  modified  system  at  low  cutoff  frequencies.  The 
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eigenvalue  analysis  showed  a large  unstable  closed  loop  pole  at  low 
cutoff  frequencies.  Thus,  at  low  frequencies  the  two  separate  analysis 
techniques  are  in  direct  conflict.  The  real  problem  is  then  to  decide 
which  technique  is  giving  correct  results. 

At  high  cutoff  frequencies  (200  -*  800  Hz) , there  is  less  differ- 
ence between  frequency  domain  and  state  space  analysis  results.  The 
eigenvalue  analysis  continues  to  show  an  unstable  closed  loop  pole 
until  cutoff  frequency  is  higher  than  725  Hz.  This  unstable  pole 
could  explain  the  erratic  nature  of  the  open  loop  stability  figures. 

It  could  also  explain  the  closed  loop  response  figures  exceeding  the 
worst  case  figures.  As  a result,  the  state  space  approach  seems  to  be 
the  correct  way  to  analyze  the  modified  system.  This  modern  control 
approach  has  become  quite  refined  and  is  particularly  suitable  to 
multi-input,  multi-output  control  system  analysis. 
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IX . Ii  ocQnr.endnti  ons 

Tnc  ADAij  control  syston  is  lo.sically  a r:i:ltiple-in;.)ut/r.’Jltiple- 
output  syste;-;.  In  the  original  ADA.i  control  syston  desiyn,  a pjrcat 
deal  of  effort  n’cnt  into  enourin,~  the  an-lo  and  translation  loo]3s 
vere  coj'.nlctely  dec.onplod  electi'onicai.ly . Ancn^  other  thinys,  this 
alloncd  frooucncy  donain  techniques  to  be  used  in  dosi(jnin(j  the 
compensation  netuor':s.  If  tt.e  vibrational  environnent  uas  well 
defined  and  T.irror-r'.otors  with  infinite  torque  and  acceleration 
capability  were  available,  tlus  technique  would  have  woir.cd.  olnce 
thesa  thj.nyr,  arc  not  avo.ilo.Llc  soiiiC  ccnprorlses  I'.ust  be  made.  The 
thrust  of  this  study  vra.s  to  iynore  the  exact  decouplin;;.:  and  decrease 
the  torque/acceleration  uci.ia’id  to  r.otor  ,/l . since  the  translational 
errors  are  not  a.c  critical  as  the  cUiq'ular  errors,  m.oi'e  translational 
error  would  be  accepted  in  order  to  decrease  this  torque/accelern,tlon 
requircr.ent . 

However,  as  soon  as  the  exact  decoupling  is  altered,  the  Units 
of  the  frequency  donain  nodcl  arc  viola.ted.  It  is  210  longer  possible 
to  assune  singlc-input/single-output  relationships,  dince  accepting 
nore  translationa.1  error  and  eliranating  exact  decoupling  rciiains  a 
viable  solution,  it  is  clean  th.at  another  modeling  ai'.d  a.na,lysis 
technique  is  necessary . 

It  is  rcconncndcd  that  modern  state  space  control  techniques  be 
applied  to  this  problem.  A state  space  model  of  ,\DA3  would  be 
capable  of  handling  multiple-inputs  and  multiple-outputs,  fince  a 
new  model  is  recommended,  it  is  further  recommended  that  completely 
nev/  feedback  and  compensation  be  designed.  In  this  way  the  resulting 
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dccipi  v.’ould  not  bo  a r.odif ication  to  the  original  control  £ycto:r.  tut 
a conplctoly  nc;  deeijn. 

To  illucti'ato  tl-iic  new  dcairpi  philosophy,  a bloc];  diagran  of  the 
basic  Airboiuio  Dynanic  Aliennent  Ayrtcir.  is  shoun  in  Fig.  16.  Ayain 
the  desij;'i  v:ill  be  perf orned  or.  a.  two  loop  control  sycten.  In  t.ho 
final  desiy;'.  the  two  basic  two  loop  systens  could  be  joined  together 
for  the  full  .ID-td  sycter. . 'lote  that  tlicre  is  no  r.irrcr  anjle  fcedbac];, 
electronic  dccouplinj,  or  cor.pennati on  netwerh  that  a.ppears  in  the 
present  control  systen . 

Tile  ope:i  loop  .state  r.atrin  equations  that  describe  t];e  syctcr.i 
are  a.s  follow.s; 


21(t)  = ['.] 


i = [=]At) 


Z = 


X(t) 


e - 


where 


[••■]  = 


0 

■l-bC7 

0 

0 


1 

0 

0 

0 


0 

0 

0 

-1407 


0 

0 

1 

0 
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.^irijular  llicplaccr.ent  of 
-lotational  Velocity  of  E 
Anpailar  Dirplacer.ent  of 
fotational  Velocity  of  ? 


0 

1^4  0? 


0 

0 


0 

0 


C 

lAc? 


'onria’icl  to  i.'otor 


'or.mand  to  fotor 


- ^ 

e(  t)  = 

Pt(  t) 

r^'cins  iii'i’oir  ol'  Bctin 

T'  ‘ 

Anjle  frror  of  fcan 

— .. 

^ J 

r -] 

— ^ 

d..,(  t) 

- 

rrceislat;  onal  fisturbanoe  of  Bear. 

1 

1 

A:i£;i;lar  fisiurbance  of  Bcai'i 

'..'hen  thf  r.  taric  ryaten  is  choun  in  a state  r.atrix  diat;ran,  it 
is  easily  identified  as  a tasic  tracivin"  prcblerr..  Tliis  state  riatrix 
cUa'^ra.i:.  is  in  7i3.  1?.  ihe  siijnificai'ice  of  identifying  this  basic 
traci'.inc  problcn  is  that  there  arc  several  syster.atic  vrays  available 
of  solving;  this  problcn.  (fef.  1:.?53). 


! 

I 
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Fip:.17’.  Tap-ic  ADA3  Trac.’-’.nr  rroblcn 


'.’ityi  the  ir.Ptfht  (gained  from  conducting  the  previous  ctudy  of 

DA'3  in  the  frequency  f^onain,  the  following  ir  an  outline  of  the 

olution  approach  for  the  rtatc  space  ADAS  tracking  problem; 

1.)  Improve  "open  loop"  mirror-motor  damping  by 
feedin'-  tho  states  (X(t))  through  a matrix 
to  the  input  u^t) . 


2.)  Ir.provc  overall  r.yrteri  clorer]  ]c;op  rcrponro  ly 
fcedlac!;  of  Iho  dlr.t’.’.rbae.cc  throiiyl.  a i''.atri:: 
to  the  input  u(t). 

lipn  the  rtatcr.  are  foci  bach  through  ncctrix  control  u(t) 

"hoconer  j 


u(t) 

nuhrtitutiny  u(t)  into  the  X(t)  equation; 

xCt)  = 


V / 4.  N 

:iv 


\\  ■ / 


1(0 


A + 


x\ 

iL\ 


Thur.,  boeor.cr-  the  n;”-'  natrtx  th.at  tlcfinon  the  clocel  loop 

rcnpor.ce  b'^twtc-n  u(t}  and  l_(t}.  A'ov-',  ’.:.lth  the  proper  choice  ef  th.c 
state  feedback  r.atri  x t;  o clor  'd  loop  eiye.nval'.ior  r.ay  be 

chanpcd  to  improve  tha:  danpiny  of  the  ryrtcn.  (."cf.  . 

"hen  the  dinturbance  (^(t))  i'-  fed  bad:  to  the  input  (2l(t))  '..'ith 
the  revised  0]’cn  loop  ryoter,  the  control  ''yctc;':  v.'ill  respond  to  correct 
any  input  dir,turb?jicc.  nio  nc?-,r  closed  loop  ryrten  c;quation  is  nov;; 

^(t) 

It  is  through  this  ratrlx  that  the  relative  inporta,r;cc  of 

the  translational  and  aiiyula.r  errors  nay  be  specified.  Cinco  the 

angular  error  is  of  jreatcr  innortance,  f K,”!  r.^y  be  designed  '.rith  a 

L J 

heavier  vfGi"htinc  on  angular-  errors.  Tiic  best  way  to  dcrien 
is  through  linear  optimal  control  tiicory'.  A coot  function  (j)  vfould 
be  defined  to  r.ako  angular  error  control  more  important  than 


II 

A + ?r 

X 

1(0  + 

‘ d 

1 

» . 

- — 
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trar.r:lational  error  control  . 


'r!ic  coot  function  (j)  in  alno  the 


place  to  define  th.c  reirhtir.f  of  the  con.trole  a arid.  u^. 

The  only  prol'lcr;  in  ir’.nlor.cnti! n^  tlic  aVovc  rccoiir.cndationo  in 
that  tuo  of  the  otatco.  and  and  Voth  din  tart -jneen  d^(t)  and 

d,(t)  arc  not  directly  nvallahlc.  A'jain,  nocern  control  theory 
provider,  a rolution  to  thin  dllen.na.  dince  two  of  the  stator  and 
arc  availahle  Un'DM,ph  ontpnt  y(-)i  '*'•  ohrerver  nay  he  used  to 
estir.ate  the  stater  X^  and  All  four  rtotes  nay  then  he  fed  tack 

thron.^h.  [^ZvJ  inprove  the  ryrten  dar.pint;.  Tlie  dirturhaneer 
(£l(t))  are  also  er-tinatod  uninr  a IXalr.an  filter,  (Tcf.  The 

Xalna-'i  filter  used  the  errors  that  arc  available  from  the  rencorr.  and 


L'o:n  the  obrcr\'er  to  estinatc  _d(  0 • finally 


d(t)  :s 


uvltiplicd  hy  into  the  control  u(t). 

of  the  total  ryrten  in  shown  in  Pi".  If. 


ooveral  .wore  detailr  need  to  he  cnaninee  before  the  control 


nyrten  shoirr  In  Pij.  If  v.’ill  work  properly,  Thc.re  incl’.idc:  pick.iny 


actual  eigenvalues  to  thr-  denired  clorcd  loop  renponso,  and 

deternination  of  weiphtiny  natricien,  on  e.^(t),  e.(t),  ii.  (t),  and 
u.,(t).  Additionally,  the  ohrerver  and  Xaln.an  filter  would  need  to  ho 
"fine  tuned"  for  optir.al  cstir.a,tcn  of  the  stater,  ^lnd  dirturhaneer . A 


conputer  proi-rar,  'would  have  to  he  designed  to  rd  nulate  this  new 
control  syrten.  It  would  then  he  through  rirr.ulation  of  th.c  rysten 
that  thc.re  details  could  he  worh.ed  out,  Cuch  a desiyn  approach  is  felt 
to  he  advantageous  hecau.ro  it  is  inherently  capable  of  dcalinp;  with 
ti;e  nulti-lnput  and  nulti -output  .system.  Tor  this  particular  rysten, 
it  docs  not  appear  beneficial  to  exactly  decouple  and  treat  as  tvfo 
sinyle-input/sinylo-output  syrten  dcsiyns. 
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Appendix  A 


Frequency  Domain  Computer  Program 


ADAS  is  modeled  in  the  frequency  domain  through  a general  purpose 
frequency  domain  program  called  "FRQPSP."  The  total  program  consists 
of  22  subroutines  and  a small  main  program.  Ihe  function  of  the  main 
program  (FRQPGM)  is  to  initialize  the  run  and  call  the  main  subroutine 
"FRQRSP."  The  calling  sequence  is  shovm  below. 


(1)  (2)  (3)  (4)  (5)  (6) 


FRQRSP 

FIRST 

FUNCTIONG 

GET 

LAST 

FiraCTIONPSD 

PRINT 

SAVE 

SOLVER 

CCEF 


SIEDEL 

SPARSE 

PIVSEL 

STEP 

FRPLOT 


COEFADS 

CON ADS 

CCEFAPT  BLOCKADAS 

CON  APT 

COEFAAS  BLOGKAPT 

CONAAS 


\ 
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^ ^ 

Subroutines;  FRQRSP,  SAVE,  GET,  SOLVER,  SPARSE,  PIVSEL,  SIEGEL, 

FRPLOT,  SIEP,  and  PRINT  belong  to  the  general  purpose  package. 

Subroutines;  COEF,  CCEFADS,  GONADS,  BLOGKADAS,  COEFAPT,  GONAPT, 

BLOCKAPT,  COEFAAS,  CCNAAS,  and  FUNGTIONG  are  used  to  describe  the 
APT,  AAS,  and  ADAS  subsystens  to  the  ALL.  All  of  these  subroutines 
were  retained  in  the  simulation  during  this  study  to  facilitate  later 
total  system  integration.  COEFADS,  GONADS,  BLOGKADAS,  and  FUNCTTONG 
are  the  only  subroutines  that  wore  altered  for  the  ADAS  control  system 
' modification.  These  subroutines  are  listed  at  the  end  of  this  appendix. 

) I ' 

; f I 

The  "FRQRSP  USER  GUIDE",  (Ref.  3)  contains  a detailed  description  of  ; 

the  rest  of  the  subroutines.  It  also  gives  the  user  the  information 
I necessary  to  adapt  FRQRSP  to  any  particular  problem. 

Output  is  in  tabular  and  graphical  form.  Frequency  is  in  both 
radians  per  second  and  Hertz.  Amplitude  is  in  magnitude  of  G(s)  and 
' decibels,  where  one  decibel  equals  20  Log^^  | KAC(G)|.  The  phase  of 

I G(s)  is  tabulated  in  degrees  as  the  imaginary  part  of  G(s) . Also, 

f ihase  is  corrected  to  account  for  phase  beyond  + 180  degrees. 
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Appendix  B 

Diagrams  of  Two  for  One  Beam  Steering  Movement 

The  understanding  of  the  two  for  one  mirror  movement  is  crucial 
in  understanding  the  ADAS  control  system.  The  need  for  the  2:1  ratio 
is  discussed  in  the  main  section  of  this  report  and  again  in  Ref.  1. 

This  appendix  attempts  to  explain  the  need  through  the  use  of  several 
diagrams.  These  diagrams  are  not  of  the  actual  layout  of  the  beam 
steering  mirrors  but  the  principle  still  applies. 

Fig.  19  shows  the  Initial  position  of  each  mirror,  the  sensor 
location,  and  the  nominal  position  of  the  optical  beam  on  target.  The 
beam  is  completely  aligned  along  its  nominal  path.  In  Fig. 20  , the 
beam  has  a pure  translational  misalignment  X at  the  sensor  face.  In 
order  to  move  the  beam  to  strike  the  center  of  the  target  both 
mirrors  must  be  moved  counterclockwise  exactly  the  same  amount.  Mirror 
#1  must  always  move  to  such  an  angle  that  the  beam  strikes  the  hinge 
(pivot)  point  of  Mirror  §2,  This  hinge  point  is  called  the  "critical 
hinge  point."  Fig.  21  depicts  the  beam  with  pure  angular  misalign- 
ment at  the  sensor  face.  To  correct  for  this  pure  angular  disturbance 
Mirror  #1  must  move  twice  as  far  as  Mirror  #2.  In  a continuous  system, 
Mirror  #1  would  also  be  required  to  move  twice  as  fast  as  Mirror  #2. 

It  is  possible  to  have  both  angular  and  translational  misalign- 
ment of  the  beam.  This  condition  is  shown  in  Fig.  22  . For  the 
combination  of  errors,  the  mirror  movement  ratio  (1:2)  will  be  some- 
where between  2:1  and  1:1.  The  actual  ratio  will  depend  on  the  relative 
magnitude  of  each  error  involved.  In  Fig.  22  , the  ratio  is  1.8:1. 
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This  principle  of  bean  steering  movement  may  be  directly  applied 
to  the  ADAS  beam  steering  mirrors.  The  only  difference  between  these 
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figures  and  the  physical  beam  steering  mirrors  is  that  the  errors 
may  also  occur  due  to  movement  of  the  target. 
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Fig. 23  . ADA3  State  Space  Hodel  Computer  Frogram 
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Fig.  25"  State  Space  Matrix  [a]  for  Modified  ADAS  T;.-o  Loop  Control  Syatc: 


F1g.27.  Comparison  of  TranKlaUonal  ADAS  SC.i  and  Filter  PSDs 


FiC.28.  ADAo  3C'.J  PCiD  and  State  Space  Filter  Approx ir.ati on 
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alter  the  closed  loop  response.  This  is  due  to  the  fact  that  the  cutoff 
frequency  of  the  filter  must  be  extremely  high  to  keep  the  modified  system 
stable.  It  is  concluded  that  the  complimentary  filter  does  not  solve  the 
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